Ab initio studies of spin-spiral waves and exchange interactions in 3c? transition metal 

atomic chains 



J. C. Tungi and G. Y. Guo^'^* 
^Department of Physics and Center for Theoretical Sciences, 
National Taiwan University, Taipei 106, Taiwan 
^Graduate Institute of Applied Physics, National Chengchi University, Taipei 116, Taiwan 

(Dated; January 12, 2013) 

The total energy of the transverse spin-spiral wave as a function of the wave vector for all 3ci 
transition metal atomic chains has been calculated within ah imtio density functional theory with 
generalized gradient approximation. It is predicted that at the equilibrium bond length, the V, 
Mn, and Fe chains have a stable spin spiral structure, whilst the magnetic ground state of the Cr, 
Co and Ni chains remains to be coUinear. Furthermore, all the exchange interaction parameters 
of the 3d transition metal chains are evaluated by using the calculated energy dispersion relations 
of the spin-spiral waves. Interestingly, it is found that the magnetic couplings in the V, Mn and 
Cr chains are frustrated (i.e., the second near neighbor exchange interaction is antiferromagnetic), 
and this leads to the formation of the stable spin-spiral structure in these chains. The spin- wave 
stiffness constant of these "id metal chains is also evaluated and is found to be smaller than its 
counterpart in bulk and monolayer systems. The upper limit (in the order of 100 Kelvins) of the 
possible magnetic phase transition temperature in these atomic chains is also estimated within the 
mean field approximation. The electronic band structure of the spin-spiral structures have also 
been calculated. It is hoped that the interesting findings here of the stable spin-spiral structure and 
frustrated magnetic interaction in the 3d transition metal chains would stimulate further theoretical 
and experimental research in this field. 

PACS numbers: 71.70.Gm, 73.21. Hb, 75.10.-b, 75.30.Ds 



I. INTRODUCTION 



Noncollinear magnetism, especially the spin-spiral 
structures, has received much attention in recent decades, 
not only for possible magnetism-based technolo gica l 
applications [ll-y| but also for fundamental physics pl-[l3|. 
In particular, it was recently reported that spin chi- 
rality in geometrically frustrated pyrochlore compounds 
could generate magnetic monopoles^ ^] and also large 
anomalous Hall effect [l0[. It was also proposed recently 
that the spin-spiral structure could be the main source 
of the magnetoelectric effect observed recently in multi- 
ferroic oxides. 0, Q The spin-spiral structure, in which 
the magnetization rotates along a certain direction in 
a bulk material, was observed two decades ago in neu- 
tron diffraction experiments on fee Fe and Fcioo^ Co^: 
alloy precipitates in Cu.[13] This experimental finding 
has since stimulated many ah initio studies of the spin- 
spiral structures in bulk magnets [l^, [13413 ■ Indeed, ah 
initio calculations [H, [13, [l3] corroborated that stable 
spin-spiral states exist in fee Fe. Furthermore, ah ini- 
tio total energy calculations for the spin-spiral structures 
also helped to formulate an explanation of the anomalous 
magnetovolume properties of the Invar alloys (the Invar 
effect) [3 • 

Noncollinear magnetism in low-dimensional systems 
has also been studied both theoretically and experi- 
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mentally in recent years. 0, [2l| - [23| For example, Mn 
monolayer on W(OOl) surface was recently investigated 
jointly by spin-polarized scanning tunneling microscopy 
and also ab initio calculations, and it was concluded 
that a spin-spiral structure along the (110) direction 
exists in this monolayer system. A stable spin-spiral 
structure with propagation vector q = (0,0,0.15)(27r/a) 
was also predicted to exist in the unsupported free- 
standing Fe(llO) monolayer with the lattice constant of 
3.16 AgU, HI]. Co/CoPt bilayers were also found to 
support noncollinear spiral structures by Brillouin light 
scattering [23"]. Interestingly, very recent ah initio calcu- 
lations show that in Mn chains on Ni(OOl), the magnetic 
structure could change from non-coUinear to collinear fer- 
rimagnetic, depending on whether the number of the Mn 
atoms is even or odd. 6] 

Stimulated by possible unusual magnetism in one- 
dimensional (1-D) systems, we have recently carried out 
systematic ah initio studies of the collinear magnetic 
properties of linear and zigzag atomic chains of all 3(i[2^. 
4d and Srifisj transition metals. Although the ideal infi- 
nite freestanding id transition metal atomic chain is un- 
stable and cannot be prepared experimentally, short sus- 
pended monostrand metal nanowires and atomic chains 
have been prepared in mechanical break junctions. [26l - [28| 
Furthermore, structurally stable Co atomic chains have 
been recently prepared on a vicinal surface of Pt(997) 
surface[29] or inside nanotubes.[30| Therefore, we have 
also performed ah initio calculations for the 3c? transition 
metal linear atomic chains on the Cu(OOl) surfacc[31, 32j 
in order to understand how the substrates would affect 
the magnetic properties of the nanowires. 
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The purpose of the present work is to study possible 
spin-spiral structures in all 3d transition metal atomic 
chains by ab initio calculation of the total energy of the 
spin-spiral state as a function of propagation wave vector 
q. Indeed, we find that the magnetic ground state in the 
V, Mn and Fe chains would be a spin-spiral state. Fur- 
thermore, we evaluate the exchange interaction parame- 
ters between the atoms and also spin-wave stiffness con- 
stants of all the atomic chains considered here from the 
calculated energy dispersion relations of the spin-spiral 
waves. The obtained exchange interaction parameters 
allow us to understand why the spin-spiral state is stable 
in the V, Mn and Fe chains but is not stable in the Cr, 
Co and Ni chains. Finally, we also estimated the upper 
limits of the magnetic phase transition temperature for 
all the atomic chains. 

This paper is organized as follows. After a brief de- 
scription of the computational details in Sec. II, we 
present all the calculated energy dispersion relations of 
the spin-spiral waves of the 3d atomic chains in Sec. III. 
These results show that a stable spin-spiral state exists 
in the V, Mn and Fe chains. Reported in Sec. IV are the 
obtained exchange interaction parameters which enable 
us to understand the stability of the obtained magnetic 
ground state in each atomic chain considered. In Sec. 
V, we present the calculated spin-wave stiffness constant 
and also the estimated magnetic phase transition tem- 
peratures for the 3d atomic chains. Finally, the band 
structures of the spin-spiral state of the V and Mn chains 
are displayed in Sec. VI, and a summary of this work is 
given in Sec. VII. 



II. THEORY AND COMPUTATIONAL 
METHOD 

In the present first principles calculations, we use the 
accurate frozen-core f ull-p otential projector augmented- 
wave (PAW) method, [33| as implemented in the Vienna 
ab initio simulation package (VASP) [U,!!^. The calcu- 
lations are based on density functional theory with the 
exchange and correlation effects being described by the 
generalized gradient approximation (GGA)[36|. A very 
large plane-wave cutoff energy of 500 eV is used. The 
shallow core 3p electrons of the 3d transition metals are 
treated as valence electrons. We adopt the standard su- 
percell approach to model an isolated atomic chain. The 
nearest wire and wire distance adopted here is 20 A. We 
start with the theoretical equilibrium bond lengths for 
coUinear magnetic states from our previous study of 3d 
TM nanowires.^24j However, in the fully unconstrained 
noncoUinear magnetic calculations [st'] for the spin spiral 
structures, we vary the bond length in order to study the 
bond length dependence of the stability of the spin spi- 
ral state. The F-centered Monkhorst-Pack scheme with 
a fc-mesh of 1 x 1 x n (n = 100) in the full Brillouin 
zone (BZ), in conjunction with the Fermi-Dirac-smearing 
method with a = 0.02 eV, is used for the BZ integration. 



TABLE I: Calculated equilibrium bond length d, ground state 
spin-spiral wave vector (g) and total energy [E{q)] [relative to 
that of the FM state {q = 0)] as well as spin magnetic moment 
(rris) at g = in the 3d transition metal chains. 







d 


q 


E{q) 




(/is / atom) 


(A) 


(27r/d) 


(meV/atom) 


V 


1.47 


2.05 


0.25 


-166.4 


Cr 


4.18 


2.32 


0.50 


-155.2 


Mn 


4.43 


2.40 


0.33 


-113.3 


Fe 


3.30 


2.25 


0.10 


-15.5 


Co 


2.18 


2.15 


0.00 




Ni 


1.14 


2.18 


0.00 





We consider the transverse spin-spiral states where all 
the spins rotate in a plane perpendicular to the spiral 
propagation vector q. The total energies of the transverse 
spin-spirals as a function of the magnitude of spin-spiral 
wave vector q are calculated self-consistently by using the 
generalized Bloch condition approach [sl, H^. To study 
the exchange interactions, we apply the frozen-magnon 
approach and obtain the exchange interaction parameters 
by a Fourier transformation of the energy spectra of the 
spin-spiral waves. 



III. STABILITY OF SPIN-SPIRAL STATES 

The calculated total energies [E{q, 9)] as a function of 
the spin-spiral propagation vector q of the 3d transition 
metal chains at several different bond lengths d are plot- 
ted in Fig. [TJ Since we consider here the transverse 
spin-spiral waves only, the angle between the chain axis 
(i.e., z-axis) and the magnetization direction = 7r/2, 
and hence we simply write E{q,6 = tt/2) = E{q). The 
spin-spiral structure at wavevector vector q = corre- 
sponds to the collinear ferromagnetic (FM) state, whilst 
the state at g = 0.5 (27r/(i) corresponds to the antiferro- 
magnetic (AF) state. Therefore, as shown in Fig. [TJ at 
q = the lowest total energy state of the Cr (Mn) chain 
occurs at 2.80 (2.60) A, but it appears at 2.32 (2.29) A 
a,t q = 0.5 {2i: /d), being in good agreement with our pre- 
vious collinear magnetic calculations 2J| . Interestingly, 
Fig. [1] shows that both the FM and AF states in the 
V, Mn and Fe chains become unstable against formation 
of a spin-spiral structure. Furthermore, the lowest total 
energy of the spin-spiral state occurs at the bond length 
that is generally different from that of the collinear mag- 
netic states. For example, the ground state of the Mn 
chain is the spin-spiral state with the equilibrium bond 
length of 2.40 A, instead of 2.60 A (the FM state) and 
2.29 A (the AF state) [H]. 

Nonetheless, there is no stable spin-spiral state in the 
Cr, Co and Ni chains. In the Cr chain, therefore, the 
AF state remains the stable state. In the Co and Ni 
chains, the FM state still has the lowest total energy (see 
Fig. [T|). In fact, we could not even obtain a spin-spiral 
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FIG. 1: (color online) Total energy E{q) versus spin-spiral 
wave vector q of the 3d transition metal chains at several 
different bond lengths d. 



solution for the Ni chain at the wave vector q being larger 
than 0.3 (27r/(i). This is because, as shown in Fig. [1] an 
increase in the number of valence electrons leads to an 
increased stabilization of the FM state whilst a decrease 
in the number of valence electrons tends to stabilize the 
AF state. This observation is further corroborated by 
the fact that in the Cr chain of bond length d = 2.32 
A, the total energy decreases steeply as the spin spiral 
wave vector increases (Fig. [T}d). For comparison, we 
notice that in previous GGA calculations pOl| . the FM 
state could not be stabilized in bulk Cr metal, whilst 
the magnetization energy of the AF state is rather small 
(-0.016 eV/atom). 

The calculated equilibrium bond length (d), ground 
state wave vector [q) and total energy [E{q)] as well as 
spin magnetic moment per unit cell (tos) at g = of all 
the 3d transition metal chains considered here are listed 
in Table H] We notice that the spin magnetic moment 
of the spin-spiral state of the V chain is significantly re- 
duced, in comparision with the FM state (see Ref. H^), 
but is close to that of the AF state. Interestingly, Ta- 
ble H] shows that the stable spin-spiral wave length A of 
the Fe chain is nearly exactly 10 bond lengths (or lat- 
tice constant), whilst, in constrast, that for the V and 
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FIG. 2: (color online) Calculated spin-wave energy spectra 
e{q) [i.e., magnon dispersion relations hhj{q)] of the 3d tran- 
sition metal atomic chains at the ground state bond length. 



Mn chains is much shorter, being 4 and 3 bond lengths, 
respectively. However, the spin-spiral wave length A can 
depend on the bond length, and this dependence is es- 
pecially pronounced for the V chain, as demonstrated 
in Fig. [1^. When the V chain is stretched to the bond 
length of 2.40 A, the spiral propagation vector q becomes 
~0.35 (27r/(i), but when it is further stretched to d — 2.55 
A, q is reduced to —0.20 {2'K/d). A similar behavior of 
the spin-spiral wave vector can be found for the Fe chain, 
as shown Fig. [T]i. 

The energy of a spin- wave excitation (i.e., the magnon 
dispersion relation) is given as the derivative of the total 
energy of the ^in-spiral state with respect to the magnon 
number [19, kll 



s{q) = huj{q) = 2fii 



' AM 



rrisQ 



[E{q)-Em (1) 



where AE{q,9) = E(q,e) - E{{),9) = E{q) - E{0) is 
the energy of a spin spiral of wave vector q relative to 
the ferromagnetic state {q — 0), AM is the decrease of 
the magnetization per site projected onto the z-axis, and 
niso is the spin magnetic moment per site at q = 0. The 
calculated magnon dispersion relations for the V, Cr, Mn, 
Fe, Co and Ni chains at the minimal energy bond length 
are plotted in Fig. [5] 

Ab initio calculations of the excitation energy of the 
spin-spiral wave along some high-symmetry lines in the 
Brillouin zone in bulk bcc Cr^20iJ, fee Mni23^^ FefH- 
m, i^, fee Co[iMl [13], and fee NifMM HO] met- 
als have been reported before, and the stable spin-spiral 
structures were found in bcc Cr and fee Mn. The sta- 
ble spin-spiral wave with a wave length of —7 lattice 
constants was also found in a freestanding bcc Fe(llO) 



TABLE II: Calculated exchange interaction parameters (Joj) 
(meV) between two j*'' near neigbors {j — 1, 2, 3, 4, 5) in the 
V, Cr, Mn, Fe, Co and Ni atomic chains. 





Joi 


Jo2 


Jos 


Jo4 


Jos 


V 


4.2 


-22.8 


-0.2 


2.6 


-2.4 


Cr 


-65.0 


17.6 


-7.2 


4.0 


-2.6 


Mn 


-78.4 


-43.2 


9.8 


-3.0 


1.2 


Fe 


158.2 


-57.8 


4.4 


2.6 


-3.2 


Co 


156.4 


13.0 


-22.8 


16.6 


-11.8 


Ni 


109.0 


6.6 


20.0 


-24.6 


7.4 



monolayer with lattice constant of 3.16 A in two previous 
ab initio studies[2ll[23[. Here we predict the existance of 
the stable spin-spiral structures in ID freestanding tran- 
sition metal (V, Mn and Fe) atomic chains for the first 
time. 



IV. EXCHANGE INTERACTIONS 

To a rather good approximation, we can map a metallic 
magnet onto an effective Heisenberg Hamiltonian with 
classical spins [Ulil 



(2) 



2J 



where is an exchange interaction parameter between 
atomic site i and site j, and tTj (crj) is the unit vector rep- 
resenting the direction of the local magnetic moment at 
site i {j). In the frozen magnon approach, the exchange 
interaction parameters Jij are related to the magnon ex- 
citation energy £(q) by a Fourier transformation 



Joj - 



-iq R 



^(q) 



(3) 



where Nq is the number of q points in the Brillouin zone 
included in the summation and 
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Here, 9 is fixed to tt/2 for all the spin-spiral states and 
J(q) is the Fourier transform of the exchange parameters. 
We therefore evaluate the exchange interactions in the 
Sd transition metal chains via Eqs. [3] and U by using the 
calculated magnon dispersion relations, as shown in Fig. 

m 

The obtained exchange interaction parameters as a 
function of the interatomic distance are plotted in Fig. 
|3] and also listed in Table HIl In the minimum energy 
bond lengths, as shown in Fig. [3^, the magnetic cou- 
pling between two first nearest neighbors in the V, Fe, 
Co and Ni chains is ferromagnetic (Jqi > 0), whilst it is 
antiferromagnetic ( Jpi < 0) in the Cr and Mn chains. In 
the Co and Ni chains, the magnetic coupling between the 



J(q) 



123456789 
Interatomic distance (bond length) 

FIG. 3: (color online) Calculated exchange interaction pa- 
rameters Joj versus interatomic distance for the 3d transition 
metal chains with different bond lengths as labeled. In (a), 
the atomic chains are in their respective ground state bond 
lengths, and in (b), the atomic chains are in either stretched 
or compressed bond lengths. 



second nearest neighbors remain ferromagnetic ( J02 > 0) 
and this explains why the ground state of these chains 
is ferromagnetic. In contrast, the magnetic coupling be- 
tween the second nearest neighbors in the V, Mn and Fe 
chains is antiferromagnetic (J02 < 0), i.e., the exchange 
Interactions in these chains would be frustrated. As a 
result, noncollinear spin-spiral states in these chains may 
become energetically more favorable than either coUinear 
ferromagnetic or antiferromagnetic state. Let us take the 
Mn chain as an example, and consider the atomic spin 
at the origin. This spin tends to couple antiferromag- 
netically both with its two nearest neighbors and also 
with its two second near neighbors because the first and 
second near neighbor exchange parameters Jqi and J02 
are negative (see Fig. [3^ and Table [TT|) . However, this 
would make its two nearest neighbor spins "frustrated" 
because they would have to couple ferromagnetically with 
one nearest neighbor on one side and antiferromagneti- 
cally with the other nearest neighbor on the opposite 
side. This frustrated magnetic coupling therefore would 
energetically favor a spin-spiral state. In fact, accord- 
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ing to the mean field theory for an 1-D classical Heisen- 
berg spin chain with the negligibly small magnetic cou- 
pling between third near neighbors and beyond (see, e.g., 
Ref. 112), in the frustrated magnetic coupling situation 
(Jo2 < 0), the system would be ferromagnetic (g = 0) if 
Joi > 4| J02I and would be antiferromagnetic {q = tt/2) if 
Joi < — 4| J02I. Table HIl shows clearly that the condition 
Joi > 4|Jo2| is fulfilled for the Co and Ni chains, giving 
rise to the ferromagnetic ground state, and that the con- 
dition Jqi < — 4| J02I is met for the Cr chain, resulting in 
the antiferromagnetic ground state. Interestingly, a sta- 
ble spin-spiral structure with the spiral propagation vec- 
tor q given by cos{qd) = —J±l^J2 would occur if J02 < 
and 4Jo2 < Joi < 4|Jo2|.[42| Using the exchange cou- 
pling parameters listed in Table HIl we would obtain the 
spiral propagation vector q — 0.24, 0.32 and 0.13 (27r/(i), 
respectively, for the V, Mn and Fe chains. These esti- 
mated q values agree very well with that obtained by the 
fully selfconsistent total energy calculations (Table 

Figure 1 shows that the stability of spin-spiral state in 
the V, Mn and Fe chains depends pronouncedly on bond 
length d. Therefore, we also calculate the exchange in- 
teraction parameters for these chains with several other 
bond lengths and the results are plotted in Fig. 3(b). The 
equilibrium bond length for the FM V and Mn chains is 
2.60 A. The equilibrium bond lengths for the AF Mn 
and Fe chains are 2.29A and 2.15A, respectively. It is 
clear from Fig. 3 that the calculated exchange interac- 
tion parameters can be sensitive to the bond length. For 
example, the magnitude of the nearest neighbor exchange 
interaction Jqi in the V atomic chain is dramatically in- 
creased from 2.1 meV to 40 meV as the bondlength is 
increased from 2.05 A to 2.60 A. A similar behavior is 
found for the Fe chain (Fig. 3). Interestingly, in con- 
trast, the magnitude of the nearest neighbor exchange 
interaction Jqi in the Mn chain is significantly reduced as 
the bond length is increased. Nevertheless, the exchange 
interaction parameters for the second near-neighbor and 
beyond are less affected by the bond length (see Fig. 3). 

Ab initio evaluation of the exchange interaction pa- 
rameters in bulk ferromagnets Fe, Co and Ni have been 
reported many times before. For example, the nearest 
neighbor exchange interaction Jgi in bcc Fe was deter- 
mined to be 39.0 meVp^ and 57.3 meVpsj. respectively. 
The nearest neighbor exchange interaction Jqi in fee Co 
was estimated to be 29.5 meVflSl. In bulk fee Ni, the 
Joi was calculated to be 5.6 meV|18{. In the freestand- 
ing square Fe and Co monolayers with the Cu(OOl) lat- 
tice constant, the nearest neighbor exchange interaction 
Joi is 92.5 and 77.6 meV 11], respectively. These indi- 
cate that in 2D systems, in general, the exchange interac- 
tions arc significantly enhanced as compared with their 
bulk counterparts, mainly because of reduced coordina- 
tion numbers. Table II show that the nearest neighbor 
exchange parameters Joi in the Fe, Co and Ni chains, are 
much larger than in their bulk counterparts. This may be 
attributed, at least partially, to the fact that the nearest 
bond length in the atomic chains is significantly shorter 



than in their bulk counterparts '245. O^'" calculated Joi 
values in the Fe and Co chains are also significantly larger 
than the corresponding Jqi values in the freestanding Fe 
and Co monolayers [nj. Interestingly, in the Fe atomic 
chains deposited on the Cu(117) surface, the effective ex- 
change interactionparameter Jf.f / was determined to be 
around 136 meV[44l|, being comparable with the corre- 
sponding Joi value of the freestanding Fe chain listed in 
Table II. 



V. SPIN- WAVE STIFFNESS AND CRITICAL 
TEMPERATURE 

The calculated energy dispersion relations of the spin- 
spiral waves e{q) for the id transition metal chains at 
the ground state bond length are displayed in Fig. O In 
the range of small g, e(g) = Dq^ , where the spin- wave 
stiffness constant D relates the spin-wave energy huj(q) 
to the wave vector q in the long wavelength limit. The 
spin-wave stiffness constant D of an atomic chain can 
be estimated by fitting an even order polynormial to the 
corresponding spin-wave spectrum shown in Fig. [51 The 
spin-wave stiffness constant D obtained in this way for 
the 3d metal chains are listed in Table IIIII A negative 
value of D means that the FM state is not stable against 
a spin-spiral wave excitation. Table IIIII shows that the 
spin-wave stiffness constant D in the V, Cr, Mn and Fe 
chains is negative. Only in the Co and Ni chains, the D 
is positive. 

In principle, one can also calculate the spin-wave stiff- 
ness constant D viafl™ 



D = 



2^iB d^E(q) _ 
rrisQ dq^ 3m 



(5) 



where Joj are the exchange interaction parameters and 
Roj = |Ro — Rj I is the distance between site and site j. 
In practice, Eq. |5l cannot be used directly to obtain reli- 
able values for the spin- wave stiffness constant, because 
the numerical uncertainties at the long distances are am- 
plified by the factor i?oj.[l8] Here we use this expression 
to understand the calculated Ds listed in Table Hill For 
example, the magnitude of the spin-stiffness constant D 
of the V chain is much larger than that of the Cr chain 
because the V chain has a much smaller spin magnetic 
moment (see Table I) and also a negative second near 
neighbor antiferromagnetic exchange parameter (see Ta- 
ble^. Furthermore, even though Jqi is positive in the 
V and Fe chains, the D is negative, because the V and 
Fe chains have J02 < and Jqi < 4|Jo2|. 

For comparision, the spin-wave stiffness constants Ds 
for the three-dimensional (3D) and two-dimensional (2D) 
Fc, Co and Ni systems from previous ab initio calcula- 
tions and experimental measurements, are also listed in 
Table |TTT1 It is clear from Table |llll that the spm-wave 
stiffness constant tends to become smaller as the dimen- 
sionality of the system gets reduced. This may be ex- 
pected (see Eq. 151) because the number of near neighbors 
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TABLE III: Calculated spin wave stiffness constant D 
(meVA^) and magnetic transition temperature Tc of the 3d 
metal chains. Also listed are the Ds and Tcs for the 3D and 
2D metal systems from previous ah initio calculations and 
experimental measurements, for comparision. 





stiffness D 




Tc (K) 






3D 


2D 


ID 


3D 


2D 


ID 


V 






-424 






94 


Cr 






-106 


31P 




414 


Mn 






-504 






274 


Fe 


250", 330' 


164" 


-78 


1414", 1043'* 


1265" 


410 


Co 


663", SIO' 


570", 427" 


616 


1645", WSS** 


1300" 


606 


Ni 


756", 555-', 




656 


397", 627"' 




458 



"Theoretical calculations (Ref. (1^). 

'Neutron-scattering measurement extrapolated to K (Ref. [45 
"Theoretical calculations (Ref. ) . 
■^Experimental measurements (Ref. [4^ ). 
"Brillouin light scattering measurement (Ref. [47||). 



■^Neutron-scattering measurement (Ref. 
'Neutron-scattering measurement (Ref. 



I). 



becomes smaller as the dimensionality of the system gets 
reduced. Among the three bulk 3d elemental ferromag- 
nets, Fe has the smallest spin-wave stiffness, and, inter- 
estingly, the stiffness of Fe becomes negative when its 
dimensionality is decreased to one (Table IIIip . 

Within the mean-field (MF) approximation, the criti- 
cal temperature {Tc) of the magnetic phase transition for 
the effective Heisenberg Hamiltonian can be estimated 
via the approximate expression (sol Isij ] 



:^(q) 



(6) 



where q is the spin-spiral wavevector and J(q) is the 
Fourier transform of the interatomic exchange parame- 
ters 







In the ferromagnetic case (q = 0) [s^l , 



T rpMF ™sO 



(7) 



(8) 



Using the calculated exchange interaction parameters Jqj 
(Table Hill and Fig. [3^), and Eqs. [6] and [3 we estimate 
the transition temperatures T^^ for the 3d transition 
metal chains, as listed in Table IIIII Encouragingly, the 
ferromagnetic transition temperatures evaluated by us- 
ing Eq. [8] for the Co and Ni chains are 622 and 444 K, 
respectively, being in good agreement with that obtained 
by using Eqs. in]and[7](Table|nD). Tablelmlindicates that 
the critical temperature for the 3d atomic linear chains 
considered varies from several tens to a few hundreds 
of Kelvins. The critical temperatures for the 3d atomic 
chains are much smaller than the corresponding ones for 
the bulk metals and their monolayers (Table HlI)) . This 



is because of the much reduced coordination numbers in 
these ID systems, although the exchange interaction pa- 
rameters in the atomic chains are generally larger than 
their counterparts in the bulks and monolayers. 

The more accurate expression for the critical tempera- 
ture within the random phase approximation (RPA) also 
exists [Hill 



1 



6HB 1 



1 



M N, ^ e(q) 



(9) 



where A is the magnetic anisotropy energy. The RPA ex- 
pression generally gives the critical temperatures for bulk 
ferromagnets Fe, Co and Ni being in better agreement 
with experiments than the MF expression [l^. Further- 
more, the critical temperatures with the MF approxima- 
tion are usually significantly higher than that from the 
RPA calculations. [ill, [l^ Using our calculated spin- wave 
dispersion relations (Fig. 2) and also the A values from 
Ref. m we obtain T^^^ = 230 K for the Co chain and 
T^P^ = 229 for the Ni Chain. Therefore, the estimated 
critical temperatures listed in Table ITTll should be consid- 
ered only as the upper limits. Finally, we note that in 1- 
D isotropic Heisenberg model with finite-range exchange 
interactions, there is no spontaneous magnetization at 
any nonzero temperature because fluctuations become 
important. [5^ Nonetheless, this discouraging conclusion 
has to be revised in the presence of a magnetic anisotropy 
and long range interactions. Indeed, ferromagnetism in 
1-D monoatomic Co metal chain on a Pt substrate has 
been recently reported. [1^ A detailed discussion on pos- 
sible finite temperature spontaneous magnetization in 1- 
D systems has been given in Ref. l53|. 



VI. ELECTRONIC BAND STRUCTURE 

In order to study how the spin-spiral structure affect 
the electronic band structure [3£] and also to help get- 
ting further insight into the spin-spiral instability at the 
microscopic level. [gllTsI the electronic band structures of 
the V and Mn chains at several spin-spiral wave vectors q 
are displayed in Fig. |3]and Fig. [SI respectively. The fer- 
romagnetic (g = 0) band structures are presented in Fig. 
[Jjand Fig. [5^. Because of the linear chain symmetry, the 
bands may be grouped into three sets, namely, the nonde- 
generate s- and dz2 -dominant bands, doubly degenerate 
(Axz^ dyz) bands, and (dj.2_j,2 ,dxy) dominant bands, see 
Figs. [3^- and [5^. The (dj.2_j,2, d^y) bands are narrow 
because the dj.2_y2, and dxy orbitals are perpendicular 
to the chain, thus forming weak 7 bonds. The d^z, and 
dy^bands, on the other hand, are more dispersive due to 
the stronger overlap of the dxz and dyz orbitals along the 
chain, which gives rise to the tt bonds. The s- and dz2- 
dominant bands are most dispersive since these orbitals 
form strong cr bonds along the chain. 

Two main changes could appear in a ferromagnetic 
band structure when a spin-spiral wave is introduced. 
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FIG. 4: (color online) Band structures of the V chain in differ- 
ent spiral propagation vector q along F to ^ direction. Figure 
(a), q = {2Ti/d), that indicates a spin-polarized ferromag- 
netic band structure whilst figure (f), q — 0.5 {2n/d), indi- 
cates an AF band structure. 



FIG. 5: (color online) Band structures of Mn in different spiral 
propagation vector q along T to Z direction. Figure (a), q = 
(27r/d), that indicates a spin-polarized ferromagnetic band 
structure whilst figure (f), q — 0.5 {2n/d), indicates an AF 
band structure. 



VII. CONCLUSIONS 



namely, the lifting of the accidental degeneracy at the 
cross-point of spin- up and spin-down bands and the " re- 
pulsion" of opposite-spin bands, [sl] These changes can 
be clearly seen in Figs. |3Jd and[SjD. For example, in Fig. 
2)3, the spin-up and spin-down d^^ and dyz bands clearly 
move away from each other. Of course, these changes due 
to the noncoUinear spin-spiral wave become more pro- 
nounced as q increases (Figs. |3]and[S|). In the V chain, 
the repulsion of opposite-spin bands appear to lower the 
spin-up dxz-dyz band considerably while in the mean time 
raise the spin-down d^z and dyz band significantly (see, 
e.g.. Figs. HJd and c). This could be the reason why the 
spin-spiral structure is energetically favored over the fer- 
romagnetic state. In the Mn chain, on the other hand, 
as the wave vector q becomes nonzero, the fully occu- 
pied spin-up dxz-dyz band near the F point is pushed 
down substantially (see Fig. [5]). This change could lead 
to a lowering of the total energy, thereby stabilizing the 
spin-spiral structure in the Mn chain. Figs. S] and [S] 
also show, in contrast, that as the spiral wave vector q 
becomes nonzero, the dispersive s- and 1^22 -dominant va- 
lence bands near the F point are pushed up in energy, and 
these shifts are more or less proportional to the modulus 
of spin spiral vector q. This upward movement of the s- 
and -dominant bands might raise the total band en- 
ergy. Therefore, the final equilibrium spiral wave vector 
q would be determined by a trade-off of these two con- 
trasting changes in the electronic band structure of the 
atomic chain. 



We have calculated the total energy of the transverse 
spin-spiral wave as a function of the wave vector for all 
Zd transition metal atomic chains within ah initio den- 
sity functional theory with generalized gradient approx- 
imation. As a result, we predict that at the equilibrium 
bond length, the V, Mn, and Fe chains have a stable 
spin spiral structure. Furthermore, all the exchange in- 
teraction parameters of the Zd transition metal chains are 
evaluated by using the calculated energy dispersion rela- 
tions of the spin-spiral waves. Interestingly, we find that 
the magnetic couplings in the V, Mn and Cr chains are 
frustrated (i.e., the second near neighbor exchange inter- 
action is antiferromagnetic) , and this leads to the forma- 
tion of the stable spin-spiral structure in these chains. 
The spin-wave stiffness constant of these 3cZ is also eval- 
uated and compared with its counterpart in bulk and 
monolayer systems. We have also estimated the upper 
limit of the possible magnetic phase transition tempera- 
ture in these atomic chains within the mean field approx- 
imation. The electronic band structure of the spin-spiral 
structures have also been calculated. We hope that our 
findings of the stable spin-spiral structure and frustrated 
magnetic interaction in the id transition metal chains 
would stimulate further theoretical and experimental re- 
search in this field. Indeed, after learning our ab initio 
results, Sandvik recently studied a spin-1/2 Heisenberg 
chain with both frustration and long-range interactions 
by exact diagonalization. [54] He found a first-order tran- 
sition between a Neel state and a valence-bond-solid with 
coexisting critical fc = 7r/2 spin correlations, [s^ 
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